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0^ . Abstract 
ON ' 

■ Linear quantum cellular automata were introduced recently as one of the models of quan- 

tum computing. A basic postulate of quantum mechanics imposes a strong constraint on any 
quantum machine: it has to be unitary, that is its time evolution operator has to be a unitary 
transformation. In this paper we give an efficient algorithm to decide if a linear quantum cellular 
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automaton is unitary. The complexity of the algorithm is 0(n r + 1 ) = 0(n 3 ) in the algebraic 
computational model if the automaton has a continuous neighborhood of size r, where n is the 
size of the input. 
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1 Introduction 

Q\ ■ The classical models of computation, such as Turing machines, random access machines, circuits, 

or cellular automata are all universal in the sense that they can simulate each other with only 
polynomial overhead. These models are based on classical physics, whereas physicists believe that 
the universe is better described by quantum mechanics. 

Feynman || [|| pointed out first that there might be a substantial gap between computational 
models based on classical physics and those based on quantum mechanics. The quantum Turing 
machine (QTM), the first model of quantum computation, was introduced by Benioff [l], If - Deutsch 
in |H described a universal simulator for QTMs with exponential overhead. Bernstein and Vazirani 
^ ', H were able to construct a universal QTM with only polynomial overhead. 

Other quantum computational models were also studied recently. Deutsch || has defined the 



model of quantum circuits, and later Yao [IS] has shown that QTMs working in polynomial time 
can be simulated by polynomial size quantum circuits. Physicists were also interested in quantum 
cellular automata: Biafore Q considered the problem of synchronization, Margolus 0] described 
space-periodic quantum cellular automata and Lloyd [12, [DJ discussed the possibility to realize 
a special type of quantum cellular automaton. Linear quantum cellular automata (LQCAs) were 
formally defined by Watrous Jig] and by Diirr, LeThanh and Santha 0. In the former paper it was 
shown that a subclass of LQCAs, partitioned linear quantum cellular automata (PLQCAs) can be 



simulated by QTMs with linear slowdown. Van Dam [17] defined space-periodic LQCAs and gave 
a universal instance of this model. 
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We should make clear at this point that most of these models are only theoretically motivated. 
Real life quantum computers as built today in laboratories are essentially small quantum circuits 
or partitioned cellular automata. 

A quantum computational device is at any moment of its computation in a superposition of 
configurations, where each configuration has an associated complex amplitude. A superposition 
is valid if it has unit norm. If the device is observed then a configuration will be chosen at 
random, where the probability of a configuration to be chosen is equal to the squared magnitude 
of its amplitude. Therefore it is essential that valid superpositions be transformed into valid 
superpositions, or equivalently, that the time evolution operator of the device preserve the norm. 
This property is called well-formedness, and thus it is a natural problem to decide if a given quantum 
machine is well-formed. In the case of QTMs and PLQCAs there exist easily checkable constraints 
on the finite local transition function of the machine which are equivalent to its well-formedness. 
Such constraints were identified respectively by Bernstein and Vazirani || and by Watrous |18| . 
In the case of LQCAs no such local constraints are known, still Diirr, LeThanh and Santha T|7j 
gave a polynomial time algorithm to decide if an LQCA is well-formed. Part of this algorithm was 
improved by H0yer [flCf l using a different approach. 

However, one of the basic postulates of quantum mechanics imposes an even stronger constraint 
than norm-preserving on the time evolution operator. It actually requires that this operator — as 
any other quantum operator — be a unitary transformation. We will call a machine which satisfies 
this constraint unitary. In [|| and |l8| it was proven that norm-preserving already implies unitarity 
in the case of QTMs and PLQCAs. It is also trivially true for machines with finite configuration 
sets, such as quantum circuits. But this is not true for LQCAs; it is quite simple to construct a 
well-formed LQCA which is not unitary]^ 

In this paper we give an efficient algorithm to decide if an LQCA is unitary. The complexity 
of our algorithm is cubic if the input LQCA has continuous neighborhood (most papers in the 
literature about classical linear cellular automata deal only with such cases). Our algorithm will 
use the procedure of @ which in quadratic time decides if the LQCA is well-formed. The present 
paper actually gives an algorithm which decides if a well-formed LQCA is also unitary. 

Well-formedness is equivalent to the orthonormality of the column vectors of the time evolution 
operator; unitarity requires orthonormality also from the row vectors. Deciding unitarity is much 
harder than deciding well-formedness. One way of seeing this is that whereas the column vectors 
have finite support, the row vectors can have an infinite number of non zero components. 

2 The computational model 

Let us fix for the paper the following notation. If u and v are vectors in some inner-product space 
over the complex or the real numbers, then {u\v ) will denote the inner product of u and v, and ||it|| 
the norm of u. If M is a matrix in such a space, then M* denotes its conjugate transpose. 

We recall here the definition of a linear quantum cellular automaton (LQCA) which is the quan- 
tum generalization of the classical one-dimensional cellular automaton. A more detailed description 
of this model can be found in and in 0]. 

An LQCA is a 4-tuple A = (E,g, N, 5). The cells of the automaton are organized in a line, 
and are indexed by the elements of Z. The finite, non-empty set E is the set of (cell-) states, and 
g G E is a distinguished quiescent state. The neighborhood N = (a\, . .. ,a r ) is a strictly increasing 

X A well-known example is the classical LCA Xor= ({0, 1}, 0, (0, 1), 5), where S(x,y) = (x + y) mod 2. This cellular 
automaton is injective for finite configurations but not surjective, thus the associated LQCA is well-formed but not 
unitary. 
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sequence of integers, for some r > 1, giving the addresses of the neighbors relative to each cell. 
This means that the neighbors of cell i are indexed by i + ai, . . . ,i + a r . An automaton is simple 
if its neighborhood is an interval of integers, that is a r = ai + r — 1. In this paper we deal only 
with simple automata, and we will only explain briefly in the conclusion how our results apply to 
the general case. 

The states of the cells are changing simultaneously at every time step according to the local 
transition function. This is the mapping 5 : Y,\ N \ — ► C E , which satisfies that for every (xi, . . . , x r ) G 
S r , there exists y G X such that [5(xi . . . x r )](y) / 0. If at some time step the neighbors of a cell 
are in states x\,... ,x r then at the next step the cell will change into state y with amplitude 
[5(x\ . . . x r )](y) which is denoted by (8(xi . . . x r )\y). The quiescent state q satisfies that 

The set of configurations is where for every configuration c, and for every integer i, the state 
of the cell indexed by i is Cj. A configuration c is finite if its support {i : ci ^ q} is finite. We 
are dealing only with LQCAs which evolve on finite configurations. We will denote the set of finite 
configurations by Ca, and from now on we use the word configuration to mean a finite configuration. 
For a configuration c, let idom(c) be the interval domain of c, which is the smallest integer interval 
containing the support of c. For the sake of definiteness, we define the empty interval as [0,-1] 
which is the interval domain of the everywhere quiescent configuration. 

The local transition function induces the time evolution operator which we write in matrix 
form Ua ■ Ca x Ca — > C, where Ua(cI, c) is the transition amplitude of changing configuration c to 
configuration d in one step. It is defined by 

U A (d,c) = Y[(S(ci + N)\di), 

i<=Z 

where (5(cj+7v) is a short notation for 5(ci +ai , . . . ,Q +ar ). This product is well-defined since c has 
finite support. 

The automaton evolves on superpositions of configurations which are elements of the Hilbert 
space £2(C A ). If at some time step the automaton is in the superposition u 6 C Ca , then at the next 
time step it will be in the superposition U A u. Therefore U A is also an operator on Ii{Ca)- A is 
well-formed if U A is norm-preserving, and we say that it is unitary if U A is a unitary transformation. 

We will work in the algebraic computational model where complex numbers take unit space, 
and arithmetic and logical operations take unit time. The description size of an automaton is 
clearly dominated by the local transition table 5. Therefore we define the size of the automaton to 
ben= |S r+1 |. 

For the rest of the paper we will fix a well-formed simple LQCA A = (£,<?, N,S). Without 
loss of generality we assume N = (0, 1, . . . , r — 1). Indeed let A' = (X, q, N', 5) be the well-formed 
simple LQCA with the general neighborhood N' = (j, . . . ,j + r — l) for some integer j. We claim 
that A is unitary if and only if A' is unitary. Let A" = (£, q, (j),5') be the shift cellular automaton, 
where 5' is the identity. A" is unitary since 5' is unitary. Moreover U A = U A »U A > which proves the 
claim. 
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2.1 Example 

The figures in this paper will illustrate our algorithm with the following LQCA: 
Qflip= ({a, &}, a, (0, 1), 5) with (5(x, y)\z) defined for all x,y, z £ {a, b} by the table: 



xy\z 


a 


b 


aa 


1 





ba 





1 


ab 


1/V2 


1/V2 


bb 


1/V2 


-1/V2 



Using the algorithm in [Q] it can be shown that Qflip is well-formed, and in this article we show 
that its evolution operator Uq flip is even unitary. This LQCA has an interesting property. For 
n > 0, let c n be the configuration which is b in all cells of index i S [—n, —1] and is a elsewhere. 
Then we have for all n > 1, C/q FL i P (c 1 , c n ) = (l/y/2) n . Thus an infinite number of configurations 
lead with non-zero amplitude to the single configuration c . 

For the all quiescent configuration c°, we have C/q FLIP (c°, c°) = 1, thus there is a unique con- 
figuration leading to it. Therefore when the LQCA Qflip runs backwards in time, every cell with 
index i < — 1 depends on cell —1. From this we conclude that there cannot be a LQCA with finite 
neighborhood whose evolution operator is exactly Uq fup . 

This makes the model of LQCA different from QTMs, since for every well-formed QTM M, 
there exists a QTM which runs M backward in time with a constant time overhead. It explains a 
bit why it seems difficult to simulate any LQCA by a QTM. 



3 The main result 

The main result of the paper is the following theorem. 

Theorem 1 There exists an algorithm which takes a simple LQCA as input, and decides in time 

3r-l 

0(n r + 1 ) = 0{n A ) if it is unitary. 

Since in [JiJ a 0(n 2 ) algorithm is given to decide if an LQCA is well-formed, we will give only 
an algorithm which decides if a well-formed LQCA is unitary. The following lemma states that we 
only have to verify that the rows of the time evolution operator are of unit norm. 

Lemma 1 Let U € <C CaxCa be a linear operator. If U is norm-preserving then its rows have norm 
at most 1. // all the rows are of unit norm then U is unitary. 

Proof Let c be a configuration, and c the superposition which has amplitude 1 for c and 
elsewhere. Then the norm of the row indexed by c in U is ||C/*c||. Since U is norm-preserving 
||C/*c|| = ||[/[/*c|| and the projection of UU*c on c has norm 

\(c\UU*c)\ = \(U*c\U*c)\ = \\U*cf. 

But the projection of UU*c on a unit vector has norm at most ||C/*c||, and therefore ||J7*c|| < 1. 

For the second part of the lemma observe that the projection of UU*c on c has norm 1. Since 
U is norm-preserving the projection on any other basis vector c' must be 0. Thus (c'|f7J7*c) is 1 if 
c = c' and otherwise, or in other words UU* = I, which concludes the proof. □ 

The outline of the proof is the following. First we give a sequence of reduction steps in section || 
to a graph theoretical problem and to another one from linear algebra. We then give an algorithm 
to solve the problem. The different steps of the algorithm are presented in sections [|, |6| and 0. The 
proof of the main theorem is then summarized in section |8l 
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4 The reduction 



The different reduction steps are illustrated in the figures [l] to ||[ 

Our problem is the following. We have to decide if all row vectors of the evolution operator 
associated to a given LQCA have unit norm, under the assumption that the operator is norm 
preserving. The naive method fails because one would have to compute the norm for an infinite 
number of rows. Moreover for every row there can be an infinite number of non-zero entries and 
every entry is defined by a product on an unbounded number of terms. The purpose of this section 
is to reduce our problem to a finite one. 

The configuration graph is the infinite directed graph Goo(V, E) defined by V = E r_1 x Z and 
E = {((xt, i), (ty, i + 1)) : x, y E E, t G E' r_2 , i £ Z}. To our knowledge this type of graph has been 



first used by Sutner and Maas [16| to show that a particular robot motion planning problem in the 
presence of moving obstacles is PsPACE-hard. It was used again by Sutner 15] to prove that the 
predecessor of every recursive configuration is also recursive. 

A non-empty sequence (possibly infinite to the left, to the right or in both directions) of vertices 
(. . . , (wi, i), . . . ) in Gqo is a path if and only if for at at most a finite number of indices i we have 
Wi ^ q r ~ 1 and there is an edge between every two immediate vertices. Note that a sequence with a 
single vertex is already a path. We denote by F, L, R and P respectively the set of paths which are 
finite, infinite to the left, infinite to the right and infinite to both directions. Figure || illustrates a 
path of P. 




Figure 1: The configuration graph of the LQCA Qflip. The bold path corresponds to the config- 
uration . . . aababbaa .... 

We say that two paths p± and pi are compatible if the last vertex of p\ and the first vertex of 
P2 exist and they are the same. In that case the composition pi P2 is the concatenation of the 
two sequences after identifying the extreme vertices. If Pi and P2 are sets of paths, then 

Pi € Pi,p 2 G P2, 
^ l ^ 2 ' pi and p 2 are compatible 

Let d be an arbitrary configuration. It induces a weight function for the edges of Goo, where 
gd((xt,i), (ty,i + 1)) = \(5(xty)\di)\ 2 . We extend the weight function to paths and to sets of 
paths. The weight of a path is the product of the respective edge weights, and the weight of a path 
set is the sum of the respective path weights. The weight of a path consisting of a single vertex is 
1, and the weight of the empty path set is 0. We denote this weighted configuration graph by G 1 ^ 
which is illustrated in figure 0. 

Although the weight of an infinite path is an infinite product, it is well defined since all but 
a finite number of edges have weight 1. The following lemma establishes a strong relationship 
between the weight of an infinite path in and the entries of the time evolution matrix. 
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d= 



aabbabaa 



Figure 2: The configuration graph of the LQCA Qflip weighted by the configuration d = 
. . . aabbabaa .... The bold edges have weight 1, the normal edges have weight 1/2 and the dotted- 
line edges have weight 0. 

Lemma 2 For any configuration d the row indexed by d in Ua has norm y/ g<i(P). 

Proof We will show that there is a bijection h between the set of configurations Ca and the set of 
infinite paths P in such that for every configuration c and d 

g d (h(c)) = \U A (d,c)\ 2 . 

Summing up over all configurations c will immediately conclude the lemma. 
Let h : Ca — > P be defined for all configuration c by 

h(c) = (. . . , (cj . . . c i+r - 2 ,i), ■■■)■ 

Then it is a bijection, and the following equalities conclude the proof. 

9d(h(c)) = J|g d ((ci...Ci +r _2 ) i),(ci +1 ...c i+r _i,i + l)) 

= nn*(^)K*)i a 
npte+*)](di) 

= \U A (d,c)\ 2 . 

□ 

With Lemma [l] we got the following reduction of our problem. 

Corollary 1 In a well-formed LQCA, for every configuration d, we have gd(P) < 1- Moreover the 
automaton is unitary if and only if for every d, gd(P) = 1- 

Let us fix an interval [j, k] and a configuration d with idom(cZ) C [j 7 k]. This interval induces a 
subset of the path sets L, F and R. For every w, w' G S 1 '" 1 we set 

Lw = {p £ L : the last vertex of p is (w,j)}, 
pj,k _ f p . first vertex of p is (w,j), 
w ' w ' 1 ^ ' last vertex of p is (w 1 , k + 1) 
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Rw' = {p £ R '■ the first vertex of p is (u/, k + 1)}. 
Since the set of infinite paths can be decomposed as 

P= |J Li®F^ k w ,®R k w „ 

w.w'GT, 7 '- 1 



we have 

= £ • o,(J#* ,) • g d (R k w ,). 

The following lemma shows that gd(L w ) and gd(Rw) are independent from j, k and d. 

Lemma 3 For any intervals [j,k], [j',k'] and configurations d,d' such that idom(d) C [j,k], 
idom(cZ') C [j 1 , k'\ and for every w G S r_1 ; we have 

9d{L j w ) = g d '{Lt) and g d {R k w ) = 9d'{Rw)- 

Proof We prove only the first equation, the proof for the second one is analogous. Let m = f — j. 
We define a bijection from l? w to L? w which preserves the weight. If p = (. . . , (w{, i), . . . , (vjj,j)) is 
a path in L w , then we define its image as p' = (. . . , (wi,i + m), . . . , (wj, j + m)). This is clearly a 
bijection and we also have g d (p) = 9d'(p') since di = q for i < j and d\ = <? for i < j'. □ 

We define the left and n<?/t£ border vectors, respectively I = {l w ) w &T, r - 1 an d r= (r w ) w£ -£r-i as 
follows: for w S S r_1 , l w = g d (Lw) and r w = g^R 1 ^), where [j,k] is an arbitrary interval and d 
an arbitrary configuration satisfying idom(<i) C [j,k]. The next lemma states that I and r are in 

R sr "\ 

Lemma 4 For all w £ T, r ~ 1 , l w and r w are finite. 

Proof Suppose there is a to such that l w = oo. (The case r w = oo is symmetric.) We will prove 
that this implies the existence of a configuration such that the associated line vector has infinite 
norm, thus contradicting by Corollary [l] the hypothesis that A is well-formed. 

Let w' be such that r w i > 0. There exists such a w 1 , since for example r q r-i > 1. Let 
xi, X2-, ■ ■ ■ , X2r-i G S such that w = x\ . . . x r -\ and w' = x r . . . X2 r -2- We set Wi = XiXi + \ . . . Xj +r _2 
for i = 1, . . . , r and Uj = . . . £i+ r -i for i = 1, . . . , r — 1. Note that w\ = w and w r = w'. 

For i = 1, . . . , r — 1 let yi £ E be such that (5(vi)\yi) ^ 0. Let j be an arbitrary integer, and set 
k = j + r — 2. We define the configuration o! to be the quiescent state outside the interval [j, k] and 
dj+i-i = V% f° r i = 1, • • • ,r — 1. Then in already the set of paths going through the vertices 
i w i>j), ■ ■ ■ i ( w ri k + 1) has infinite weight. Since each path has non-negative weight, P has also 
infinite weight. □ 

The first part of our algorithm will be the computation of the border vectors I and r. For the 
second part, we reduce now our problem to a question in linear algebra. 

For every a E S, let M a S Mp r xSr be the linear operator whose matrix is defined for all 

w, w' £ S r_1 as 

\(5(xty)\a}\ 2 if w = xt,w' = ty for some x,y G S and i G S r_2 , 
otherwise. 



M a {w',w) 



We extend this definition to finite sequences over E. If e denotes the empty word, then M e is the 
identity operator. Let s > 1 be an integer, and b = b\ . . . b s be an element of E s . We define 

M b = M ba ■ ■ ■ M bl . 
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Lemma 5 Let d be a configuration with idom(d) = [j,k]. Then 



g d (P) = (M d] ... dk l\r). 

Proof We have 

9d(P) = E 9d(Li) ■ g d (F^ k w/ ) ■ g d (R k w ,) 
= 9d{F 3 w k w ,)-r w , 

w,w' 

= y~] l w ■ M d ._ dk (w f , w) ■ r w > 

w,w' 

= E ( ^2 lw ' M dj...d k {w',w) J • r w , 

w' \ W / 

= ^2(M dj ... dk f)(w')-r w , 

w' 

= {M dj ... dk l\r). 

□ 

Since for every b £ £*, there exists a configuration d whose non-quiescent part is b, Corollary |l] and 
Lemma [| imply the following reduction. 

Corollary 2 A well-formed LQCA is unitary if and only if for every b £ S* , we have (Mbl\r) = 1. 

We also have the following property, which simplifies our next reduction step. 
Lemma 6 For every well-formed LQCA A we have (l\r) = 1. 

Proof Let d be the all quiescent configuration. Then in Ua the column indexed by d has the entry 
1 at row d and elsewhere. Since the column vectors of Ua are pairwise orthogonal, Lemma |l| 
implies that the row indexed by d has only zero entries besides column d. Therefore this row has 
norm 1 , and the claim follows from Lemma |2| and [|. □ 

Let m = |S r_1 |. The border vectors can be seen as elements of W 71 , and the elements of the set 
M = {M a : a G £} can also be seen as linear transformations in R m . Let us fix a few notations for 
the inner product space R m . Let S C R m a finite set of vectors. The linear subspace and the affine 
subspace generated by S, denoted here respectively by (S) and [S] are defined as 

(S) = {Aisi + . . . + A t s t | i > 0; si, . . . ,s t e S; Ai, . . . , At G R}, 
[S] = {AiSi + ... + A t a t |t>0;5i,... ,s t E S;Xi, . . . ,X t G R; A x + . . . + A t = 1}. 

B is said to be a basis of (S) (respectively of [S]) if (B) = (S) (respectively [B] = [S]) and has 
minimal cardinality for this property. 

Let u £ R m be a vector and T C ]]j mxm be- a finite family of linear transformations. We set 
S + u = {v + u : v £ S}, and F(S) = {f(v) : v £ S, f G J 7 }. Let denote the linear subspace 
whose normal vector is u, that is = {v : (v\u) = 0}. 

We define by induction on i, for i > 0, the sets ^(S). Let ^{S) = S, and JP+^S 1 ) = 
F{S) U F{F(S)). We say that S is closed for v under F, if UZo^iM) Q S. 
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From Lemma ||, H r + I is the set of vectors which have unit inner product with r, that is 

H f +l = {u : (u\r) = 1}. 

It is an affine subspace since H?+ I = [Hf + 1]- Clearly for every 5 £ S*, (Mbl\f) = 1 if and only if 
Hff + I is closed for I under M . Therefore by Corollary Q our reduction steps lead to the following 
theorem: 

Theorem 2 A well formed LQCA is unitary if and only if Hp + I is closed for I under M. 

This characterization is illustrated in figure ^. 

2 



1 


r 


1 = MJ 


= M b l 







Z 




-1 


.... 








-1 





1 2 


3 



Figure 3: For the LQCA Qflip, the border vectors / and r, and the affine subspace + I. In this 
example I is a fix-point for the operators M a and Mb, which shows that t7Q FLIP is unitary. 



5 Computing the border vectors 

In this section we will give an algorithm for computing the border vectors. By symmetry, it will be 
sufficient to give it only for the left vector. The main tool in the computation will be the weighted 
border graph. Its underlying graph can be seen as a slight modification of the finite version of 
the configuration graph. This graph was also used in |JJ for checking that all the columns of Ua 
had unit norms. However, there the weights were defined as the norms of the transition state 
superpositions, whereas here they will be the squared magnitudes of the amplitude of the quiescent 
state in those superpositions. 

The (left) border graph is the finite, directed, weighted graph Gi = (V,E,g). The vertex set is 
V = S r_1 and the edge set is 

E = {(xt,ty) : x,y G S,t e Y, r ~ 2 } 
The weight function is defined as 

g((xt,ty)) = \(5(xty)\q)\ 2 . 

A path in Gi is a finite, non empty sequence of at least two vertices such that there is an edge 
between every two consecutive vertices. Observe that a single vertex alone here does not form a 
path. As usual, the weight of a path is the product of the edge weights, and the weight of a set of 
paths is the sum of the individual path weights. The weight of the empty path set is 0. 

For every w £ E r— , we define P w as the set of paths in Gi whose first vertex is q r ~ 1 , whose 
second vertex is different from q r ~ 1 and whose last vertex is w. 
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Lemma 7 For every w £ S r , we /icwe 

Lw \ g(P w ) + l zfw = q r -K 
Proof Let d be a configuration with interval domain [j, k], and let 

p 9 = (...,((r l ,o»-.((r l ,i)). 

We set L' w = L? w — {p q }- We will give a weight preserving bijection from L' w to P w which maps p to 
p'. Let p = (. . . , (wi,i), . . . , (vjj,j)) be an element of L' w , where Wj = w. Let h < j be the greatest 
integer such that for every % < h, we have = g r_1 . Then we set p' = (q r ~ 1 , w^+i, • • • , Wj). This 
is clearly an injective mapping, and it is also surjective since by the choice of h, w^+i ^ q r ~ l . 
It is also weight preserving since the edges in p until the vertex (wh,h) have all weight 1. Since 
9d(Pq) = 1 the lemma follows. □ 



3(r-l) 

Theorem 3 There exists an algorithm which computes the border vectors in time 0(n r + 1 ). 

Proof According to Lemma |7| it is sufficient to compute g(P w ) for w £ S r_1 . The main difficulty in 
this computation is that the paths of P w are defined by a constraint which forces the second vertex 
to be different from q r ~ 1 . The solution we propose codes this constraint directly in the graph, 
which we will augment by one vertex for this purpose. Then we compute the total path weight 
from i to j for all vertices To do this, we will adapt a standard algorithm which constructs the 
regular expression associated to a finite state automaton. 

Let G\ = (V, E', g 1 ) where V = V U {sq r ~ 2 } for a letter s S, 

E' = EU{(sq r - 2 ,q r - 2 y):ye?;\{q}}, 

and g'(e) = g(e) for all edges e £ E and g'{(sq r ~ 2 ,q r ~ 2 y)) = 5'((g r ~ 1 ,q r ~ 2 y))- This graph is 
illustrated in figure |||. For every w E S r_1 , let P' w be the set of all paths in G\ from sq r ~ 2 to w. 
Clearly there is a weight-preserving bijection between P^ and P w . 




Figure 4: The graphs G[ (left-hand) and G' r (right-hand), associated to the LQCA Qflip. From 
these graphs we can compute I = and f = („). 

The border vectors have only finite components, nevertheless for their computation we have 
to extend the non-negative real numbers with oo. Let R* be this set. We define the following 
computation rules with respect to oo: 



oo + c = c + oo = oo, 
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for every c£M*, 



oo-c = c- oo = oo-oo = oo, 

for every real number c > 0, 

oo • = oo • = 0, 

and oo° = 1. We also define c* for every cGl* as Y2e?=o ° e > that ^ s 

1/(1 -c) if0<c<l, 
oo otherwise. 

Let {vi,V2, • • • , be an arbitrary enumeration of the vertices of G\. For 1 < i, j < \V'\ and for 
< k < \V'\, we define the path sets P k (i,j) as the set of paths which start in Vi, end in Vj, and all 
the other vertices in the path have index less or equal to k. Let Wk(i, j) denote g(Pk(i, j))- Then 
we claim that W k (i,j) satisfies the following recursion for 1 < i,j < \V'\, and 1 < k < \V'\ : 

w ( . ~ _ f 9'{{vi,Vj)) if (vi,Vj) e E', 
01 ,J> \ otherwise, 

W k (i,j) = W k ^ 1 {iJ) + W k - 1 (i,k)-(W k - 1 (k,k))*'W k - 1 {k,j). 

We prove our claim by induction on k. In Po(i,j) the only path is the edge between Vi and Vj 
if this edge exists. 

Assume that this equation is true for k — 1. We note that for every path of P k (i,j) there exists 
a unique integer e such that vertex v k appears exactly e times the path. Thus we can write 

oo 

Pk{i,j) = Pk-i{i,j) U I) Ph-i(i, k) (8) P fc -i(fc, k) ® ■ ■ ■ ® P k -i(k, k) ®P k -i(k,j), 

e =i v r. ' 

e— 1 

where the unions are disjoint and (8) is the path composition operator defined in section ||. By 
induction hypothesis we have 

oo 

W k (i,j) = W k ^(i,j) + (w k -i{i,k) ■ (W fc _i(M)) e ' • WVi(M) , 

e'=0 

which concludes the induction. 

This proves the correction of the following algorithm: Let m = \V'\ = |S| r_1 . Initialize Wq. 
For k = 1, ... ,m compute W k using W k -\. Finally output the border vector I defined by l(w) = 
W m {sq r ~ 2 ,w) for w ^ q r ~ l and /(g r_1 ) = W m {sq r ~ 2 , q r ~ 1 ) + 1. Proceed in similar fashion for f. 

3(r — 1) 

The complexity of the algorithm is 0(|S| 3 ^' -1 ^) = 0(n r + 1 ). □ 



6 Closed affine subspace 

In this section we will give a polynomial algorithm for the following problem. 
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Closed Affine Subspace 



Input: Two vectors l,f G M m such that (l\r) = 1 and a set of linear transformations 
M = {M a : a G £} in M m , where m = |E| r-1 . 

Question: Is H? + l closed for I under M, i.e. for all b G S* do we have G i/f + i? 

We set t = |S|. For the simplicity of notation, let = Z and let Ei = M l ({l\). Since 

H = [H], we have ^ C H if and only if [Ei] C #. Therefore we have to decide if [U^o^l ^ F - 
Dimension arguments imply the existence of a fixpoint, a set £j, such that [Ej] = [U^o^]- 
Moreover we need only to keep track of a basis of [Ei], that is a set Bi of linearly independent 
vectors, with [Bi] = [Ei]. 

-, 3r-2 

Theorem 4 There exists an algorithm which decides if H is closed fori under M. in time 0(n r + 1 ). 
Proof We claim that this is realized by the following algorithm: 

Bo := {1} 
i := 1 

while [Bi] + [BiUMiBi)] 

B i+ i := a basis of [Bi U M{Bi)] 

i :=i + l 
B :=B t 

if B C H accept 
else reject 

At every iteration dim([5j]) increases, and therefore the algorithm terminates in at most m — 1 
iterations. We prove that it is correct. 

We show [Bi] = [Ei] for every i > by induction. The statement holds by definition for i = 0. 
Suppose [Ei] = [Bi] for some i. Since M contains only linear operators, we have for any set S, 
[M(S)] = [M([S])]. Therefore [M(Ei)] = [M{B{)] which implies [E i+1 ] = [B i+1 ]. 

Since B is a fixpoint, that is [B] = [B U M(B)], this implies that [B] = [USo^l and the 
correctness follows. 

We now turn to the analysis of the complexity. We will build inductively the basis so that for 
all i, Bi C Bi + \. At the i-th iteration, to build Bi + ±, initially we set = Bi. Then we compute 
every vector in A4(Bi\Bi-i) and add it to Bi + \ if it is not in the affine subspace generated by 
Bi + \. At the end of the algorithm, these steps were applied to all vectors Mb for M G M and 
b G B, thus at most \M\ ■ \B\ = 0{tm) times. Computing Mb takes 0(m 2 ) with standard matrix 

multiplication, and checking affine independence takes also 0(m 2 ) with the algorithm described in 

i 

the next section. Thus the overall complexity is 0{tm 6 ). The theorem follows since t = n r +^ and 

r-l 

m = n r + 1 . □ 



7 Maintaining a basis 

In this section we give a dynamic algorithm for the following problem. We want to maintain a basis 
B of a d-dimensional linear subspace in R m , such that the following requests for a given vector 
u G M. m can be treated efficiently: 

Membership query: "Is u G (B)?" 
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Add to basis: Replace B by B U {u}. 

We can define the problem for affine subspaces as well. Fortunately the latter can easily be 
reduced to the former: Let / : W 71 — ► M m+1 be the function which maps a vector u to u' with u\ = Uj 
for i = l,... ,m and u' m+l = 1. Then every vector v satisfies v G [B] if and only if f(v) G (f(B)). 

A solution to this problem requires a tricky data structure to encode the basis. The naive way 
would be to represent B by a matrix such that its column vectors are exactly those of B and to 
apply the Gaussian elimination algorithm (see for example to check whether u G (B). This 
would require 0(md 2 ) time steps. 

Transforming the matrix in an upper triangular form is the bottleneck of this approach. The 
representation we choose for B will improve this complexity. 

Theorem 5 There is a dynamic algorithm for maintaining a basis of a d- dimensional sub space in 
M. m which treats each request in time 0(m(m — d)). 

Proof We represent a non empty basis B by the couple (T, B), where T G R mxm is an orthogonal 
matrix and (T(B)) = M. d x {0} m ~ d . The empty basis is represented by (1, 0), where I is the identity 
matrix. 

Since T((B}) = (T(B)), u G (B) if and only if Tu € R d x {0} m_d . Thus verifying u G (B) 
is reduced to checking if the last m — d components of Tu are all which can be done in time 
0(m(m — d)). 

Suppose u G" {B). We will show that there is a orthogonal matrix M affecting only components 
from d + 1 to m which satisfies {MTu)d+\ ^ and (MTu)i = for i = d + 2, . . . , m. Thus 
(MT(B)) = (T(B)) and (MT(BU{u})) = R d+1 x {o}™"^ 1 ). Therefore {MT,BU {u}) represents 
BU{«). 

We define M as the composition of two operators M\ and M2 we describe now. By hypothesis 
u G" (B), therefore there exists an index k G {d + 1, . . . , m} such that {Tu)k 7^ 0. Define M\ to be 
the permutation matrix which exchanges k and d + 1. 

Let = M\Tu. Note that 7^ 0. Then for an arbitrary vector v we define (M2v)i = Vi 
for i = 1, . . . , d + 1 and (M2v)i = Vi — Vd+iu^/u^^ for i = d + 2, . . . , m. Clearly Mi and M2 are 
orthogonal linear operators, and since 

M 2 uf = (u[,u 2 , ... , u' d+1 , 0, ... ,0) 

M satisfies the required property. We can compute M 2 M\T in time 0(m(m — d)), which concludes 
the proof. □ 

8 Putting all together 

We are now able to prove Theorem ffl, that is to give an algorithm to decide if a given LQCA is 
unitary. By Theorem || to solve this problem we have to compute the associated border vectors 
and decide the corresponding Closed Affine Subspace problem. According to Theorem a the 

3(r- 1) 

border vectors can be computed in time 0(n r+1 ), and due to Theorem B the last problem can be 

3r-2 

solved in time 0(n r+1 ), which concludes the proof. 
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9 Conclusion 



A not necessarily simple LQCA can be transformed into a simple one with the same time evolution 
operator. Let the original neighborhood be N = (oi, ... ,a r ). The size of the new neighborhood 
will be s = a r — a\ + 1. If we define the expansion factor of an LQCA as e = (s + l)/(r + 1) then 

3r-l „ 

the algorithm works in the general case in time 0(n r+1 ) = 0(n ). 

In the case of space-periodic configurations van Dam (T^j has shown that LQCAs can be effi- 
ciently simulated by QTMs. Watrous [18] gave an equivalent result for partitioned LQCAs. This 
question remains still open for the model of this paper. 
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